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Abstract 

For classical lattice systems with finite (Ising) spins, we show that the 
implementation of momentum-space renormalization at the level of Hamilto- 
nians runs into the same type of difficulties as found for real-space transfor- 
mations: Renormalized Hamiltonians are ill-defined in certain regions of the 
phase diagram. 

PACS: 64.60Ak, 05.50+q, 02.50Cv 



1 Introduction 

Despite the great success of renormalizat ion-group (RG) ideas, both for computa- 
tions and as a heuristic guide, many aspects of the theory still lack rigorous mathe- 
matical justification. The filling of this gap is more than just of academic interest. 
It has been repeatedly pointed out (e.g. || p. 82], || footnote in page 38], [[13], 
p. 268]) that the method is not a black-box type of technique; its succesful appli- 
cation requires some understanding of the underlying physics or one may be led 
to incorrect conclusions. Studies on the foundations of real-space transformations 



15| , pL6|, [L9| , f29| , |30|| suggest that a similar remark applies to the underlying mathemat- 
ics. Indeed, these studies show that in various occasions renormalized Hamiltonians 
are ill-defined. The finite-volume probabilities of the renormalized system exhibit a 
long-range dependence on boundary spins that is incompatible with the existence of 
a Hamiltonian, at least one defined in the usual (summable) sense. Such a "pathol- 
ogy" is usually referred to as non-Gibbsianness . This phenomenon, which appears 
after a single application of the transformation, was first detected in the vicinity of 
first-order phase transitions, but was later discoverd in other regions of phase dia- 
grams, including at high magnetic fields |SD|, and at high temperatures ^E, f23 



It follows that the design of the renormalization transformation is crucial for the 
very existence of a renormalization flow in a suitable space. 

Nevertheless, the lack of similar studies for momentum-space transformations 
left open the possibility that they could be free of this pathological behavior. That 
is, the question remained as to whether such transformations, possibly with a soft 



cutoff, would generally lead to a honest-to-God renormalized Hamiltonian |I6| [|. In 
this note we present a simple example showing that this is in general not the case, 
as already suspected by Griffiths |14| |. There is no essential difference between real- 



space and momentum-space transformations at least in case the spins are bounded. 
The same mechanism — the existence of a phase transition in the system of original 
spins, constrained by a particular block-spin configuration — causes similar problems 
with the definition of renormalized Hamiltonians. For an earlier suggestion that 
momentum-space maps are not all that different from real-space maps, see also []17| . 



As we remark at the end of Section [|, these problems can be interpreted as 
a manifestation of the well-known "large-field problem". It might be hoped that 
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the considerable experience accumulated in the treatment of this type of problems, 
could be of help to control the non-Gibbsianness "pathologies" . 



2 Momentum transformations 

We consider Ising spins a% = — 1, +1 on a lattice, x £ Z d . For each finite periodic 
cube V = [-N, N] d in ~K d we define the Fourier-transformed variables 



a 



V \ _ „—ikx 

x&V 



where kx := k\X\ + • • • + k^Xd, and each hi belongs to the Brillouin zone: Bn 
{-7t, — tt[1 - 1/(2N + 1)], . . . ,tt[1 - 1/(2JV + l)],7r}. The inverse of Q is 

1 



(2A^ + l) c 



E^e*", (2.2) 



for x G V. 

A momentum-space transformation is defined in two steps: 



1st. A cutoff is applied to the variables a¥: 

a'Y := f(k) a\ . (2.3) 



The volume-independent cutoff function / is designed so as to keep only mo- 
menta smaller than a certain threshold fco- 

2nd. Momenta are rescaled by the factor k Q so as to return to a Brillouin zone in 
[-vr,7r] rf : 

a'Y := f(k%/n)al kQ/7r . (2.4) 

In addition, the renormalized variables a'Y are usually rescaled in applications. We 
will not do this, as we shall not apply the transformation more than once. 

In Wilson's original approach |31], and references therein], the cutoff function 
f(k) was chosen simply as the step function 

It was quickly realized, however, that such a sharp cutoff leads to unwanted long- 
range terms in the renormalized Hamiltonian (see e.g. pl| , p. 153], |21|, Appendix 
2]). To avoid such terms one usually takes smooth momentum-cutoffs, that is, 
functions / which go to zero in a sufficiently differentiable fashion. Such functions 
are obtained, for instance, via a convolution 

f(k) = £ 5 A (k - f) X k (k) (2.6) 
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with a smooth delta-like function 5& peaked at k = of width A. 

Rigorously speaking, one is interested in the limit V — > Z d of this procedure. 
To make sense of this limit we return to real space, where the prescription 
translates into the relation 



a$ = J2f V (Lx'-y)a g , x'eV/L, (2.7) 

y&V 

where f v is the V^-dependent (inverse) discrete Fourier transform of / and 

L := i . (2.8) 

The volume V is assumed to be a disjoint union of cubes of side L (i.e. N is a 
multiple of L). As V — > Z d , the function f v (x) tends to 

/(*) := T^p f_J(k)e-^dk. (2.9) 

A sharp cutoff in momentum space gives rise to a non-summable function /, i.e. 
Sxgz d 1/0*01 = 00 • [The inverse transform of Xfc 1S proportional to the function 
Ilf=i sin(k Xi) / (koXi)]. Summability is restored if / is smooth enough (for example 
once differentiable). In such a case, expression ( |2.7|) remains valid in the thermody- 
namic limit: 

a' sl = Yl f(Lx f -y)ay, (2.10) 

and the renormalized spins remain bounded in this limit. They may take a large 
number of values, but within some finite interval. 

Expression fl2.10| ) shows that a cutoff in momentum space, even a smooth one 



like ( p.6|) , leads to non-local averages in real space, i.e. to functions / extending to 
infinity. This is the distinctive feature with respect to the real-space transformations 



analyzed, for instance, in ||30|| . Nevertheless, it is expected that "the physics behind 
integration over fluctuations having wavenumbers [\ki\ > k ] is the same as the 
physics behind the formation of blocks of spins having volume [L d ] in real space" 
p3| , Section 4.2]. To ensure this, the momentum-space cutoff should lead to an 
almost local average in real space. That is, the function / should decay rapidly 
outside of a region of size not much larger than L. We see that if / has a Fourier 
transform of the type ( |2.6| ), the contribution of spins outside the volume of size L 
is of order ln(A//c ). We conclude that the cutoff function / must approach zero in 
a "gradual" manner, that is, with A of the order of ko in ( |2.6|) . 

A soft momentum- cutoff is a function / that is smooth and gradual in the above 
sense. 
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3 The phenomenon of non-Gibbsianness 

A state (probability measure or distribution) is called Gibbsian if it can be written in 
terms of Boltzmann-Gibbs weights for an "acceptable Hamiltonian [which] . . . must 
satisfy the additional requirement of locality . . . [that is,] a quantity that is additive 



over distant lattice sites" fl31| , p. 145]. In other words, Hamiltonians must be such 
that the energy of disjoint volumes is additive except for boundary terms whose 
contribution is small in comparison with the volumes. For classical lattice systems, 
the appropriate requirement is that the flipping of one spin lead to a finite energy 
change whatever the configuration of the remaining spins is. This intuition has been 
precisely formalized in the form of a summability condition that we now present. 
Hamiltonians are, in general, sums of many-body terms: 

H(a) = ^$ B M , (3.1) 

B 

where each $^ is a function only of the spins in the finite set B C Z d , i.e. of the 
variables Ob = {c"x}fes- F° r Ising spins, these functions $b are usually written in 
the form Jb^b] the general expression ( |3.1| ) is more suitable for spins larger than 
1/2, where one would need powers of ag, and also for some particular spin- 1/2 
interactions P7| |. The finiteness of the single-flip energy change translates into the 



following summability requirement: 

sup J2 \\®b\\ < oo , (3.2) 

* BBx 

where ||$b|| := sup^ |$s(cr B )|. 

If this summability holds, then we can consider Hamiltonians with arbitrary 
boundary conditions. For each cube A in 2 d , let 

H A {a\u) = J2 ®b((ctaw)b) , (3.3) 
where a\u is the configuration 

((TA ^ = { Z if 1 1 Z d \ A • (3 - 4) 



For finite-range interactions, the dependence of the Hamiltonian (|3.3j ) on the bound- 
ary condition u is only through spins at sites not farther from A than the range of 
the interaction. More generally, even if the Hamiltonian involves terms with arbi- 
trarily long range, the summability condition ( |3.2| ) implies that the dependence on 
the boundary spins decays with their distance to the volume A. This property, which 
is called quasilocality , is central for our argument. Lets us state it precisely. We take 
a sequence of cubes U D A with larger and larger radius and fix a configuration u 
in the intermediate ("buffer") region U\A (see Figure ffl). It is not hard to see that 



September 2, 1998 



5 





— ► 
U -^7L d 



Figure 1: Test for quasilocality. As U tends to 7L d while keeping fixed the intermedi- 
ate configuration u, the energy inside A should asymptotically become independent 
of the configuration r\ or fj outside U 



the summability condition implies (in fact, it is equivalent to) the following fact: As 
U tends to the whole ~K d the Hamiltonian becomes independent of what happens 
outside U : 



sup 

vv 



H k (a\uurj) - H K (a\uJurj) 







(3.5) 



for all o"a and all boundary conditions u, for all cubes A. Hamiltonians H\ satisfying 
( |3.5| ) are said quasilocal. 

The condition ( |3 .5| ) is appropriate for discrete finite spins. The renormalized 
spins considered in the next section, however, can take an infinite number of values. 
In such a situation the quasilocality condition must allow small variations in the 
spin configurations a and uo [|TJ . 

The Hamiltonians (|3.3|) are used to construct the Boltzmann-Gibbs weights 



Pa((t\u) :-- 



exp[— pH\(a\u) 

Z A (u) 



(3.6) 



[Z\(u) is the obvious normalization factor], which do inherit the quasilocality prop- 
erty @: 



sup 

IV 



p A (cr|o; [ /?7) - pA(a\uuv) 







v- 



(3.7) 



for all cubes A, configurations <j\ and boundary conditions uj. For spins taking 
infinitely many values, the condition should be slightly modified as described in 
note 0. 

A momentum-space renormalization amounts to passing to renormalized Boltzmann- 
Gibbs weights 



i 



Z a ,(lu>) 



exp[-PH A (a\u)] 

a A weC(cr' A/ w') 



(3.8) 
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where C{a'^,uj') denotes the set of (original) spin configurations compatible with the 
"block-spin" configuration u' k ,uj\ that is the set of rj such that 

E/(J*-flw={£ i«eZ«\A'. < 3 - 9 > 

yez d k x x 

The symbol " " m ( p.8|) reminds us that the operation involved is not a stan- 
dard sum because we are "summing" over uncountably many configurations. The 
operation is, therefore, a sum combined with a suitable limit procedure or, in 
mathematical terms, an integral with respect to the countable product measure 
Ilx e z d [(l/2) Ylt]-]- (The reader interested in the rigorous construction of (|3.8j) is 
referred to the discussion in |3y, pp. 987-90].) 

Most of the studies based on renormalization ideas assume that renormalized 
weights fl3.9| ) can also be written as Boltzmann-Gibbs weights for an acceptable — in 
Wilson's and Kogut's sense — renormalized Hamiltonian. That is, it is assumed that 
the family of identities 

exph3'ff;,(<7V)] := "]T" exp[-PH A (*\u)] , (3.10) 

gives rise to Hamiltonians H' that can be written in the form ( |3.3| ) for a suitable 
interaction that satisfies the summability condition (|3.2| ). The purpose of this paper 
is to show that this assumption can be false at least at low enough temperatures. 
This is the phenomenon of non-Gibbsianness referred to in the title of this section. 

This lack of Gibbsianness will be proven in the next section by showing that 
the renormalized weights fail to satisfy the quasilocality property ( |3.7|) , or, more 
precisely, its generalized version for spins with infinitely many values M. We shall 
show that there exists one configuration oj' for which the Boltzmann-Gibbs weight 
p(a'\(jj') violate the quasilocality condition for some volume A' (formed by only 
two sites!) and some configuration a'. We emphasize that the existence of one 
such a configuration to' is enough to disprove the existence of an acceptable (i.e. 
summable) renormalized Hamiltonian. The summability property ( |3.2j ) implies the 
quasilocality ( |3.7|) (or its generalization for spins with infinitely many vaules) for 
all u' . Nevertheless, it is natural to wonder how relevant the phenomenon is from 
the physical point of view, specially if the violation of quasilocality involves very 
atypical configurations u'. We shall comment on this point in Section f|. 

The lack of quasilocality can be interpreted as exhibiting some sort of "action 
at a distance" : Infinitely far away spin flips produce a sizeable change close to the 
origin, even when the intermediate renormalized spins are frozen in the configuration 
a/. This is in contrast with the usual behavior in equilibrium statistical mechanical 
(Gibbsian behavior) where changes at infinity can propagate only through fluctua- 
tions of intermediate spins. It is not hard to imagine the explanation: The fixing of 
a renormalized configuration still leaves some fluctuation possible in the system of 
original spins compatible with it. These fluctuations act as "hidden degrees of free- 
dom" that in some instances can bring information from infinity. This happens when 
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the constrained system of original spins develops long-range correlations, i.e. when 
it undergoes a phase transition. The argument of next section consists, precisely, in 
showing conditions and discussing examples under which such phase transitions do 
take place. 



4 Non-Gibbsianness due to momentum transfor- 
mations 

We consider the nearest-neighbor ferromagnetic Ising model in C = Z 2 

H = - £ oa^ (4.1) 

(S,y) y 

at low temperatures, that is large j3. It has been shown that the low-temperature 
states for this model under a (local) block-average transformation with even block- 
sizes are mapped onto non-Gibbsian states ]30|, Theorem 4.6]. A very simple example 
of this phenomenon for 1 by 2 blocks was presented in |p6|| . We shall now prove a 



similar result for a momentum transformation of the type introduced above. 

Let us first sketch some intuition behind our argument. In real space, a momen- 
tum transformation looks approximately like an average: 



x' 



£ f(Lx'-y)a l7 , (4.2) 



where / is the Fourier (anti)transform (|2.9| ) of the cutoff function f(k). Expressions 
like ( |4.2| ) have been studied for example in [[11], [l!J (see also |21], Appendix 2] for a 



stochastic version). Even when (O) involves a sum over all spins y of the lattice, one 
would expect that each a' s , is essentially determined only by the spins inside the 
block of side L centered at Lx' . Therefore, the mechanism causing non-Gibbsianness 
for average transformations PUL Section 4.3.5] should apply to the present case with 
minor adaptations. 

We will take for our example the identity in one direction and in the other 
direction the soft cutoff function: 

f( k ) = { cos 2 (k) for|/c| < f, 
K ' \ otherwise . 1 ' ' 

This function integrates out half of the momenta degrees of freedom in this 
direction, which corresponds to taking a (not strictly local) average over blocks of 
size 1 by 2, centered at sites with even coordinates in the direction in which we 
renormalize. 

Its Fourier transform is easily computable. Indeed, /(0) = |, f(2) = /(— 2) = |, 
and for all other n 

fin) = sin(n-) x — . (4.4) 

Jy J 7T V i> (ra - 2) n (ra + 2) y ' 
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(In particular f(n) = for all even n ^ 0, ±2). 

The initial (and crucial) part of the argument consists in exhibiting a transformed 
configuration uj' such that the corresponding constrained system of original spins has 
a phase transition at zero temperature. The configuration in question is u' 3 = for 
all x E Z 2 . The corresponding original configurations must, therefore, satisfy the 
constraint 

E /(0<4Jn+» = ( 4 -5) 
I 

for each n in the direction under consideration. We claim that the only four ground- 
states are the 4-periodic configuration (strip state) 

+ + -- + + -- + + (4.6) 

and its translates over distances 1,2 or 3 (while in the other direction they are of 
course translation invariant). It is immediate to check that these configurations are 
compatible with the constraint. Moreover, it is not difficult to check that under the 
constraint (|4.5| ) they are groundstates. 

Indeed, assume that we would have a row of 3 identical spins, say plus, next to 
each other, then we claim that the above constraint could not be satisfied. If the 
middle plus would be on an even site, say at zero, this would require that 

E/(0 > 2/(1) (=^-) (4.7) 

|J|>3 3X7F 

which a simple calculation shows to be impossible. If, on the other hand, the middle 
plus would be on an odd site, the interval of 3 plus spins would need to have a 
minus spin both to its left and to its right (because of the argument above), and the 
constraint centered at either the left or the right site could only be fulfilled if 

E /(0 > /(0) (= \) (4-8) 

|Z|>3 4 

which again is impossible. This shows that the constrained system has multiple, 
namely four, ground states. 

The remaining part of the argument follows closely the presentation in |3(J Sec- 
tion 4.3.5]. There are three additional steps: 

(1) Existence of a phase transition at nonzero temperature for the constrained sys- 
tem. This follows from a well-known theory (Pirogov- Sinai theory |24], Chapter 2], 
PD| , Appendix B]; note that as remarked in ]T3L the theory also applies to systems 
with constraints.). There is one extremal phase associated to each of the four ground 
states. 

(2) Selection of the phases of the constrained system via block-spin boundary con- 
ditions. This can be done by choosing a profile r\' such that if it is imposed in 
a sufficiently large (but finite) volume, the configuration deep inside this volume 
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has to be close to the prescribed ground state. This is straightforward, though a 
little cumbersome to write. The idea is as follows: Pick first the ground-state of 
the constrained system that corresponds to the phase one wants to select. Then 
pick a (very) large volume A, and outside it set the original spins equal to all "+" 
or all"—". Now compute the resulting block-spin configuration according to ( |4.2| ). 
It will have the property that, for regions U sufficiently inside A, the block spins 
will be very close to zero. The profile we talked about corresponds, then, to the 
block spins in some ring around A and in A \ U. One sees that by playing with the 
parity of the boundary of A and changing "+" to "— " outside A, one can select 
alternatively the four phases of the constrained system [fj]. 

(3) "Unfixing" of the spins close to the origin. This is an uncomplicated and 
inessential step. See the discussion in |30|, Section 4.2].) 



We see that the argument is insensitive to the presence of a magnetic field (be- 
cause the constrained system is asked to have small magnetization), thus we are 
proving non-Gibbsianness for low temperatures but arbitrary magnetic field. 

Moreover, we remark that in case the blocks are large, having small magnetisa- 
tion in a block at low temperatures represents a large fluctuation from the typical 
behavior, in which the magnetisation is either positive or negative of order 0(L d ). 
Renormalized effective interactions are known not to be adequate to describe such 



large values of the fluctuation field [12]. This large (fluctuation) field region can in 
good cases be treated by polymer expansions. 



5 Comments and conclusions 

The present example of non-Gibbsianness as a consequence of momentum-space 
transformations confirms the suspicion of Griffiths [14[] that " no peculiarities of this 



sort have been found. . . , which may merely reflect the fact that no one has looked 
for them!". Nevertheless, one should not draw too radical conclusions from this 
occurrence. On the practical side, the main implication of non-Gibbsianness is 
that one has to be very careful in designing renormalization group transformations. 
This is in complete agreement with what the founders and various practitioners of 
Renormalization-Group methods have been saying all along. 

Indeed, already Wilson and Kogut in their classic review emphasized 

Otherwise, [that is non-perturbatively] , the locality of [the renormal- 
ized interactions] is a non-trivial problem, which will not be discussed 
further p. 145]. 

And more explicitly, M. E. Fisher in his "Renormalization Group Desiderata", 
listed the conditions needed for a succesful renormalization scheme in Hamiltonian 
space: 
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A Renormalization Group for a space of Hamiltonians should satisfy 
the following: A) Existence in the thermodynamic limit,. . . C) Spatial 
locality . . . , one should be able to identify the same regions of space and 
associated local variables before and after the transformation || Section 
5.4.2] 

Our example adds to the numerous instances showing that perversely or sloppily 
designed transformations can lead people into trouble. As N. Goldenfeld points 



out in his book Lectures on Phase Transitions and the Renormalization Group |L3 
p. 268]: "It is dangerous to proceed without thinking about the physics" . The moral 
is, then, that renormalization transformations must be carefully crafted and case- 
tailored. Already Wilson, as quoted in |E| p. 492], warned: "One can not write a 
renormalization cookbook" . 

On the foundational side, examples like the present one confirm the view ex- 
pressed by Benfatto and Gallavotti M in the opening sentence of their book, Renor- 
malisation Group: "The notion of Renormalisation Group is not well-defined". It is 
clear that the mathematical formalization of the method requires much more than 
a naive approach in terms of Hamiltonians and flows of coupling constants. We 
mention two directions of work which can potentially lead to a better mathematical 
understanding of the renormalization group framework. 

On the one hand, the connection with the large-field problem suggests to com- 
bine renormalization group ideas with geometrical expansions — cluster or polymer 
expansions — to circumvent the ill-definedness of the renormalized Hamiltonian. 
These expansions have indeed been succesfully applied in the rigorous control of 



renormalization group transformations of unbounded-spin systems |TT] , \L% and ref- 
erences therein]. A related approach, for bounded-spin systems, resorted to the 
renormalization of Peierls-like contours [|K| . 



On the other hand, the more recent program started by Dobrushin [^|, ^|, ||, [7| 
studies non-Gibbsianness with techniques borrowed from the treatment of Griffiths 
singularities. He proposes a more general class of allowable Hamiltonians, leading to 
the notion of "weak Gibbsianness" as the right framework for a unified treatment. 

We think our result illustrates and clarifies to some extent the reason why find- 
ing a good renormalization group scheme is such a non-trivial task, not only for 
strictly local but also for only approximately local transformations. We produced 
an example in the low-temperature regime, but the fact that the mechanisms of 
non-Gibbsianness are so similar for real-space and momentum-space transforma- 
tions, leads us to the conjecture that, as in real space, also in momentum-space one 
cannot trust that in general the critical region is free of problems. 
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